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ABSTRACT 

The maximum likelihood method is used for the 
first time to restore parameters of electron- 
photon cascades registered on X-ray films. 

The method allows to carry out a structural 
analysis of the ft -quanta family darkening 
spots independent of the ^ -quanta overlapping 
degree and to obtain maximum admissible accu- 
racies in estimating the energies of the ^ - 
quanta composing a family. The parameter esti- 
mation accuracy weakly depends on the value of 
the parameters themselves and exceeds by an 
order the values obtained by integral methods. 


In our previous papers /I, 2/ we have proposed for the 
first time to use the maximum likelihood method for the 
solution of the traditional inverse problem of X-ray emul- 
sion experiments. The essence of the method consisted in 
the reconstruction of source parameters (the energy E, 
zenith and azimuth angles 0 and , the spot center co- 
ordinates X, Y) of the darkening spot on the X-ray emul- 
sion film according to the matrix of darkening densities 
{Plj } measured by a densitometer in the cells (i, j ) t ■ 
= 1, » • • • , Mi , j =1, Ni . The scanning cell dimen- 

sions AX , Ay and the number of cells hi-nj = N are 
fixed. It has been shown in /I, 2/ that the maximum likeli- 
hood methqd gives consistent, efficient and nonshifted 
estimates for the parameters X,Y in the wide range 

of scanning cell dimensions A \ » AV ■ 10,20, . .50 ju km. The 
accuracy of parameter estimation slightly depends on the 
value of the parameters themselves and exceeds by an order 
the values obtained by the momenta method and by integral 
methods. 

The results of /I, 2/ are obtained for the darkening 
spots caused by a single ^ -quantum (electron) with the 
energy E in the range 2 ♦ 100 TeV and zenith angle 
in the interval 1 ♦ 0.5* 

In /3/ the results of solution of the inverse problem 
for the general case - the number of the darkening spot 
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sources (x -quanta) is greater than 1 with the strongly- 
marked structure of the spot overlapping, are given. 

Remaining within the framework of notations accepted 
in /I ,2/, one may present the likelihood function for the 
realized density matrix {Pti } at the unknown but fixed 
vector of parameters £*, X 1# ... x*. ... Yrt, 6, V) 

in the form 


Hi Hj 


<t>({hj}la)*n n Mmi/ a)i 

1 i* 1 j=i 


( 1 ) 


where ruj * Pa • S is the number of electrons in the 
cell (i,j ); S = AX- Ay is the cell area; N i , Mj is the 
corresponding number of scanning cells along the X and Y 
axes; is the probability of H-ii electrons having 

been registered in the cell ( j ) at the given if 


J 

max 


vector 

of parameters; l ■ 1 , if Uij < m min ; £ ■ 0, if « n.^ 

i = 1 , if >a„ ax . The boundary values of the electron 

number are defined by the operation properties of the film 
and are chosen by us as in /1-3/ P m a* * n. max / S ■ 0.5jukm‘ a 
P mi = ^min / *3 88 0.002jMkm' a . The distribution function_ 

^'" n/n ) is assumed to be Poisson with the mean value n = 


a 


n. , 


a 4> _ 


X -quanta 


... , where n.^ is the number of 

(electrons) ^ a m (E m ,Xm, y m , © , y> ) is the average number 
of electrons in the cell (t, 3 ) generated by the source 
with the energy E m , whose trajectory crosses the film at 
the point X m , at the zenith and azimuth angles 0 and 9* 
respectively; n<k_ is the average number of "background" 
electrons ( = n<t >/S = 0.04ju km ). 

The determination of the q vector components is per- 
formed numerically by minimizing the likelihood function 
negative logarithm: 


•&.* min ) a> *rnin tn ^ 


( 2 ) 


The inverse problem may also be solved in the case 
when the number of darkening spot sources is unknown. 
In this case one should minimize the functional 


s "H* { mLn - a 0 ,n 0 + an,... ,u', (3) 

where n. 0 and a' determine the boundaries of search for 
the darkening source number. 

The inverse problem was solved by the maximum likeli- 
hood method basing on the density matrix {j\ j \ generated 
by the Monte-Carlo method at the given vector of parameters 
{Pij} « The average number of equivalent electrons in the 
cell (i, 3 ) was calculated according to the cascade theory 
axial approximation. Fluctuations were assumed to be Pois- 
son-distributed. The BESM-6 computer main memory allowed 
one to analyze the rasters with the cell number Hi ■ Nj $7000 
and the summary number of unknown parameters 3(hv + 1 ) $ 99. 

Consider the possibilities of the elaborated method at 
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1 ) n 1 , 2) n* » 2, that is of interest for calibrating 
the emulsion technique on JT°— 2* decays; 3) > 2 for 

processing ^ -families; 4) ly * ?, i.e. when the number of 
the darkening spot sources is the component of parameter 
vector and is unknown. 

1 . n Y = i 

The simulated rasters parameter reconstruction was 
carried out for discrete series of x -quanta energy (E 0 » 

= 2,5,10,20,50,100 TeV) A X m A y = 10,20jukm cell dimen- 
sions and various thicknesses of absorbers (d). 

The relative accuracy of estimates of E,0, made a 
few per cent. The accuracy of the definition of the spot 
center coordinates is 2-5ju km. 

2 . = 2 

For the darkening spot pair the calculations were per- 
formed at various energies E i » 2 • E a (2 ,4> 5, 1 0,20 TeV; 
and distances between the spot centers Z) » (0.1+100jm km) 
at the scanning step of 10,20,50^ km. 

The results may be formulated as follows: 

a) the accuracy of determining the zenith Q and azimuth H 
angles differs insignificantly from the single spot results 
and makes <6 CO s© ^ 5 ♦ 7 %; 

b) the accuracy of determining the spot center coordinates 
is <b x % <6u a x/2 » the distances between centers are^siAX 

c) the accuracy of determining the energy for the overlap- 
ping spots is not worse than 10% at the summary energy re- 
construction accuracy of 5-7%* 

3 • n* ^ 2 

For the overlapping families of X -quanta the inverse 
problem was solved for the energies distributed uniformly 
in the 2-20 TeV interval. The number of $ -quanta was 
varied in the region 2-20 over the area 0.25 mm • 

The accuracies of determining the parameter vector 
components at 3s 10 were independent of the ix y number 

and on the average equalled the corresponding values at 
n_ = 2. 

5 The accuracy of determining the zenith angle cosine 
made 4 + 5 % and also was independent of n-^ • 

4. ? 

The rw value varied in the interval 1 5 due to the 

complexity of solving the problem at the unknown number of 
the darkening spot sources. 

No noticeable variation of the estimation accuracy was 
observed. Note also that at the known values of rt y s 10 
the efficiency of solving the inverse problem equalled 
90 ♦ 100 %. At the unknown ru v the efficiency decreased 
to 80%, and the computer counting time grew 5-7 times. 
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